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Abstract 

We show that there are Hilbert spaces constructed from the Hausdorff 
measures 7i B on the real line M with < s < 1 which admit multiresolution 
wavelets. For the case of the middle-third Cantor set C C [0, 1], the 
Hilbert space is a separable subspace of L 2 (R, {dx) s ) where s — log 3 (2). 
While we develop the general theory of multi-resolutions in fractal Hilbert 
spaces, the emphasis is on the case of scale 3 which covers the traditional 
Cantor set C. Introducing 

^{x) = V2 Xc (3x-l) 

and 

V'zW = Xc( 3x ) - Xc( Sx - 2 ) 
we first describe the subspace in L 2 (R, (dx) s ) which has the following 
family as an orthonormal basis (ONB): 

^• fc (x) = 2*^(3^ -fc), 

where i = 1, 2,j, k 6 Z. 

Since the afhne iteration systems of Cantor type arise from a certain 
algorithm in R d which leaves gaps at each step, our wavelet bases are in 
a sense gap-filling constructions. 
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1 Introduction 

The paper has three interrelated themes: (1) construction of wavelet bases in 
separable Hilbert spaces built on affine fractals and Hausdorff measure; (2) ap- 
proximation of the corresponding wavelet scaling functions, using the cascading 
approximation algorithm; and (3) an associated spectral theoretic analysis of a 
transfer operator, often called the Ruelle operator. 

There are surprises when our results are compared to what is known for the 
traditional multiresolution approach for L 2 (M. d ), and even when compared to 
known results for special classes of affine fractals. 

Some comments on (1) - (3): Due to earlier work by Jorgensen, Pedersen, 
[TB) and Strichartz et al [23J , it is known that a subclass of the affine fractals 
admits Fourier duality. Affine fractals arise from the specification of an expan- 
sive matrix, and a finite set of translations. The fractal X itself then arises from 
this data and an iteration 'in the small' of the corresponding affine maps. Let 
L = L(X) be the associated iteration 'in the large'. We say that (X, L) is a 
Fourier duality, if an orthonormal basis on X may be built from the frequencies 
in L. While it is known that, if X is the middle third Cantor set, then there is 
no L which makes a duality pair, we show that nonetheless, every affine fractal 
admits an orthonormal wavelet basis. In our discussion of wavelets, we start 
with the middle third Cantor set; and we then pass on to the general affine 
fractals. 

As for the approximation issues in (2), we know that for L 2 (R d ), there is a 
rich family of wavelet filters which yield cascade approximation. This family of 
filters is much more restricted for the fractals: Our results for the affine fractals 
even offer a certain dichotomy ( Theorem 16. 211 : If the cascades do not converge 
in the Hilbert space, then the terms in the cascading approximation sequence 
are typically orthogonal, and thus very far from being convergent. 

Our analysis of (1) - (2) is based on spectral theory of the associated transfer 
operator, and we show in the second half of the paper (starting with Section |3J 
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how this spectral theory differs in the three cases, the standard L 2 (R d )- wavelets, 
and the special duality fractals versus the general class of affine fractals. 

Our proofs depend on ideas from geometric measure theory, and from earlier 
papers on harmonic analysis of affine fractals. While some of this material is 
in the literature, it isn't available precisely in the form we need it here. In any 
case, it is difficult for readers to locate without first having a brief overview. 
So we include a minimum amount of facts from the literature for the benefit of 
readers. We hope thereby to bridge the diverse fields, fractals, Hilbert space, 
wavelets, approximation, and harmonic analysis. 

We develop the theory of multircsolutions in the context of Hausdorff mea- 
sure of fractional dimension between and 1. While our fractal wavelet theory 
has points of similarity that it shares with the standard case of Lebesgue mea- 
sure on the line, there are also sharp contrasts. These are stated in our main 
result, a dichotomy theorem. The first section is the case of the middle-third 
Cantor set. This is followed by a review of the essentials on Hausdorff mea- 
sure. The remaining sections of the paper cover multiresolutions in the general 
context of affine iterated function systems. 

It is well known that the Hilbert spaces L 2 (M) has a rich family of orthonor- 
mal bases of the following form: 

where ip is a single function € L 2 (R), with 




and the integration refers to the usual Lebesgue measure on R. Take for example 

1>{x) = xi(2x) - XiQx - 1) (f.l) 

where / = [0, 1] is the unit interval. 
Clearly / satisfies 

2/ = /U(/ + l). 
The Cantor subset Ccl satisfies 

3C = CU(C + 2) (1.2) 

and its indicator function ip c := \c satisfies 

x 

fc(^) = fc( x ) + ( fc( x -' 2 )- (1-3) 

Since both constructions, the first one for the Lebesgue measure, and the 
second one for the Hausdorff version (dx) s , arise from scaling and subdivision, 
it seems reasonable to expect multiresolution wavelets also in Hilbert spaces 
constructed on the scaled Hausdorff measures H s which are basic for the kind 
of iterated function systems which give Cantor constructions built on scaling 
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and translations by lattices. We show this to be the case, but there are still 
striking differences between the two settings, and we spell out some of them after 
first developing the theory in the case of the middle-third Cantor construction. 

While there are other wavelet approaches to fractals in the literature, for 
example an d EH, there is in fact no overlap with this work, since the 

previous papers deal with wavelets on the fractal itself, while the present paper 
deals with wavelets on an enlarged fractal (actually a fractal measure), allowing 
a structure closer to a standard multiresolution analysis (MRA). 

The practical applications are to fractals arising in physics and in symbolic 
dynamical systems from theoretical computer science, see e.g., 0] ^7], |24|. 
There is already a considerable body of work on harmonic analysis on fractals, 
see for example [201) E3> G^L EH* an d EH- Much of it is based on subdivision 
techniques, and algorithms which use cascade constructions, but so far we have 
not seen direct wavelet algorithms and wavelet analysis for fractals. 

In section 2, we recall some facts about Hausdorff measure H s , Hausdorff 
dimension, and Hausdorff distance. They will be needed in the Hilbert space we 
build on H s . It is a natural separable subspace of the full 7Y S -Hilbert space, and 
it is built up from the algebra of Z-translations (additive), and TV-adic scaling 
(multiplicative), where N is fixed. We then turn to the cascade approximation 
for the scaling function ip defined by the usual 1/N subdivision. We prove a 
theorem for the case < s < 1 which stands in sharp contrast to the traditional 
and more familiar case s = 1 of Daubechies et. a!.; i.e., the case of the Hilbert 
space L 2 (R) based on Lebesgue measure dx on R: The scaling equation is then 

ip(x) = VN^2a k tp(Nx- k) (1.4) 

fcez 

with the masking coefficients a k satisfying the usual two axioms 

^2 Ofc = Vn, and ^ atak+Ni = <5£,o, (1.5) 
fcez fcez 

Motivated by the expression on the right hand side in i|1.4f> . we define the wavelet 
subdivision operator M by 

(Mf)(x) := VNj2akf(Nx-k), f 6 i 2 (R); (1.6) 

fcez 

and note that its properties depend on the specifications in (11.511 . 
Simple conditions are known for when the limit 

lim M n xi =<p (1.7) 

n — >oo 

exists in L 2 (R), see chapter 5. Then tp (when it exists) solves (|1.4|) . and 
there is an easy formula for building functions ipi, . . . , V'jv-i m -^ 2 (^0 from ip 
such that 

(N*i/> i (N k x-£)\l<i<N,k,£ez\ (1.8) 
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is an orthonormal basis (ONB) in L 2 (R). If N = 2, a formula for ^ is 

V>(x) = V2^(-l) fc a 1 _ fc (p(2x-/c). (1.9) 

fcez 

In general when N > 2, the functions tp 1: . . ., 4>n-i ma y result from the solution 
to a simple matrix completion problem; see ^H], [5] and |B] for details. In 
the case of Hausdorff measure Ti s (0 < s < 1, s depending on the scaling 
number N), the analogous matrix completion is still fairly simple. A main 
question (non-trivial) is now that of solving the analogue of 11. 4|) . but in the 
T^-Hilbert space. The biggest differences concern the changes in (|1.5|l and 
l|1.7fl when < s < 1. It turns out in the fractal cases that there are then 
many fewer admissible solutions than those suggested by analogy with ljl.5(l . 
We summarize the situation in Sections 4-6, where our main result takes the 
form of a dichotomy theorem; the solutions to the 7i s -convergence question 
are isolated within a larger family of masking coefficients analogous to i|1.5|) . 
There is further a new set of orthogonality conditions entering the analysis 
when < s < 1, which are not present in the more familiar case of s = 1. 

We interpret the wavelet filters as functions m on the torus T. If the scaling 
number N is given, following J3|> we introduce the wavelet-transfer operator 

(R mo f) (z) := 1 KMI 2 f{w), for / € C(T), and z e T. (1.10) 

W N =Z 

Our dichotomy for wavelets will be explained in terms of the spectral prop- 
erties of R mo j a ls° called the Ruelle operator. For simplicity, we introduce the 
normalization i? mo (l) = 1, where 1 denotes the constant function 1 on T. A 
probability measure v on T is said to be invariant if vR mo = v. Equivalcntly, 

j R ma (f)du = J fdu, for all / G C (T) , 

T T 

or 

u(R mo (f))=u(f). (1.11) 

We introduce a notion of (mo , iV)-cycles for (|l.l()(l which explains the solutions 
v £ Mi (T) to l|l.ll|l . In our setting, the dichotomy boils down to two cases for 

(i) v = Si (the Dirac mass at z — 1), or 

(ii) some v is a singular measure on T with full support. 

In the first case (i), the Hilbert space is L 2 (R); i.e., that of the standard wavelets; 
and in the second case (ii), the Hilbert space is built from Hausdorff measure 
H s , < s < 1. 

We begin the discussion with N = 3 and s — log 3 (2). 

Definition 1.1. We define 1Z to be the set of all real numbers that have a base 
3 expansion containing only finitely many ones. It is an inflated version of C: 
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TZ:= {Efel-m a ^" k \meZ,a k e {0, 1, 2} for all k eZ,a k ^l for all but 
finitely many indices k}. 

For the fractal cases, < s < 1, the factor y/N in equations H1.4(l - (|1.6(l will 
be different, see details below. Similarly the factor N? in (|1.8|l changes: With 
scaling number N and with p subdivisions, the ONB corresponding to p. 8(1 in 
L 2 (TZ,H S ), s = \og N (p) is {pz tp i (N k x — £)}. However, the geometric properties 
of the cascade approximation to the scaling function change completely as the 
Hausdorff dimension moves from s = 1 to the open interval < s < 1 . This will 
be spelled out in the last four sections of the paper. 

2 The Hausdorff Measure 

Returning to the middle-third Cantor set C = C3; i.e., N = 3 and p = 2, here 
are some elementary properties of 1Z: 

Proposition 2.1. ( The middle-third Cantor set.) The set 1Z has the following 
properties: 

(i) Invariance under triadic translation: 

k 

K+ — =K, (fc,neZ). 

on 5 v 3 ' 

(ii) Invariance under dilation by 3: 

3 n K = K, (neZ). 

(Hi) The middle-third Cantor set C is contained in 1Z and moreover it covers 
1Z by translations and dilations: 

!Z=\J\j3- n (C + k). (2.1) 

Proof, (i) Any triadic number t = with fco, no G Z, fco > 0, has a finite 
expansion in base 3: 

m 

t= J2 ^ 3 ~ fc ' t*e{o,i,2}. 

k=— m 

Take x G 1Z. Then x has a finite number of ones in its expansion so the same 
affirmation will be true for x + t. 

(ii) is clear: multiplication by 3 means a shift in the base 3 expansion. 

(iii) Since 

C= |^a fc 3- fc I a fc G {0,2}|, (2.2) 
it is obvious that CcR. 
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The inclusion "d" follows from (i) and (ii). Now take 
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x e TZ, x = y~] a k 3 k , 

fc=— m 

only finitely many are equal to 1. 

Let £ be the last index for which a; = 1, and take k : = max{n ,£o}- 
Then 

x e 3- k " \ C+ J2 a k3~ k+ko ) 

\ k——m / 

and this shows that the other inclusion is also true. | 

Remark 2.2. TZ has Lebesgue measure 0. Indeed, this follows from proposition 
2 (iii), because C has Lebesgue measure 0, and so do all the sets 3~ n (C + k) 
with n, k G Z. 

Even though some of the properties that we need for the Hausdorff measure, 
for fractals, and for iterated function systems (IFS) are known, we found that the 
material is wildly scattered throughout the literature; and to increase readability 
we have included some highpoints from these areas. This should also help bring 
out the contrast between the traditional MRA-analysis, and the present more 
stochastic approach. 

Next we define a measure on TZ. It is the restriction of the Hausdorff measure 
H s with s = log 3 (2) to TZ. 

We recall some background on the Hausdorff measures from : 

For a subset E of K., s > 0, and S > 0, define 

£ 117,1" I U t/j D E, \Ui\ <5> 

i=l i=l ) 

where \U\ = sup{|x — y\ | x, y G U} (the diameter of U). 
It is known that 7i| is an outer measure on K. 
Define 

H S (E) = limW(E) = supWf(£). (2.3) 

<5>0 

Then verify that Tt s is also an outer measure. By the Caratheodory construction 
P, if we restrict Ti. s to the cr-field of W s -measurable sets, we get a measure called 
the Hausdorff measure. 

Proposition 2.3. : 

(i) All Borel sets are measurable. 

(ii) [Inner regularity] Any Ti s -measurable set of finite Ti s -measure contains a 
F a -set of equal Ti s -measure. 

(iii) If E C TZ then there is a Gg-set containing E and of the same Ti s -measure. 

(iv) For s < 1 and G open, Ti s (G) = oo (the measure is not regular from above). 
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(v) Translation invariance: For any Ti s -measurable set E , and any t € TZ, E + t 
is H s -measurable, and 

n s (E)=n s (E + t). (2.4) 

(vi) For any Ti s -measurable set E and any c > 0, cE is Ti s -measurable and 

H s (cE) = c s H s (E). (2.5) 

Consider now Tt s with s = log 3 2 restricted to the 7i s -measurable subsets of 
TZ. We will keep the notation Ti s for the restriction. 

Proposition 2.4. : 

(i) If E C TZ is an Ti s -measurable set and t = is a triadic number, then 
E + t C R is H s -measurable and 

U S {E) = H s (E + t). 

(ii) If E C TZ is an Ti s -measurable set then 3E C 1Z is Ti s -measurable, and 

H S (3E) = 2H S {E). (2.6) 

(Hi) If f £ L 1 (1Z,TL S ) then the function on TZ, x — > /(§ ) is also in L X {JZ, H s ) 
and 

J f(x)dH s (x) = \J f(^)dH s (x). 
n ii 
(iv) Ti s (C) — \, where C is the middle-third Cantor set. 

Proof, (i) and (ii) are direct consequences of propositions 12 . 31 and 12 .41 

(iii) follows from (ii) when / is a characteristic function of an 7i s -measurable 
set. Then, for arbitrary /, the formula can be obtained by approximations by 
simple functions. 

(iv) See [Tj theorem 1.14. | 

Remark 2.5. The measure H s on TZ is still non-regular from above. All open 
sets in TZ still have infinite measure. 

To see this, we show that Tt s (I) = oo, where / = (0, 1) fl 7£. 
Indeed / D C, so H S {I) > 1. 

Also, observe that 3/ = / U (/ + 1) U (/ + 2) disjoint union (we neglect some 
points that have 7i s -measure 0.). 

Therefore, with propositions [O] (i) and (ii), we obtain 

2H S (I) = m s (I) (2.7) 

so Ti s (I) is either or oo. cannot be from the previous argument, hence it 
must be oo. 

By scalings and translations, it can be proved that TL s ((a,b) fl TZ) = oo for 
any interval (a, b) . 

Since all open subsets of TZ have measure oo it follows that no non-zero 
continuous function on K is integrable! (Just take / _1 ((a, b)) for some interval 
that doesn't contain and intersects the range.) 
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Definition 2.6. We denote by H the Hilbert space 

H := L 2 (K,H S ). 
The linear operator T on H defined by 

Tf(x) = f(x-l), (fEH,xEll), (2.8) 
is called the translation operator. The linear operator U on H defined by 

Uf(x) = ±=f (|) , (fEH,xETZ) (2.9) 

is called the dilation operator. 

From Proposition 12.41 and using some simple computations, we obtain the 
following proposition: 

Proposition 2.7. : 

(i) T and U are unitary operators. 

(a) utu- 1 = T 3 . 

Denote by ip — Xci the characteristic function of the Cantor set C. We 
prove that tp satisfies all the properties of a scaling vector. 

Proposition 2.8. The following hold: 

(i) [The scaling equation] U^p — -^(p + T 2 ip). 

(ii) [Orthogonality of the translates] (T k ip | ip} = (5/., (k € Z). 
(Hi) [Cyclicity] spMi{U n T k ip \ n E Z, k £ Z} = H . 

Proof (i) XJip = ^x 3C , T 2 (p = Xc+2 , but 3C = C(J(C + 2), so (i) follows. 

(ii) T k (p = Xc+k so T k p and <p are disjointly supported for k ^ 0. For k = 0, 
I ¥>) = J K Xc^ s = 1) by proposition EHl(iv). 

(iii) First take E C R measurable and with H S (E) < oo. We want to approxi- 
mate Xe by linear combinations of functions of the form U n T k ip. 

Note also that 

U n T k <p = X3HC+k) , (n,keZ). (2.10) 

With proposition l2.4l fiii). and partitioning E if necessary, we may assume that E 
is contained in a set of the form 3" (C + fco)- Applying dilations and translations 
we may further assume that E C C | 

Define 

V: = jr,.,, „ : - 3- n C + J2a k 3- k \ a k e {0, 2} n > 1 1 . (2.11) 

This family V is a Vitali class for E; i.e., for each x G E and each 5 > 0, there 
is a U € V with x Ell, and < < <5. 
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Indeed, we see that for all n > 1: 

Cn,a„,...,ai = C 

ai,...,a„G{0,2} 

Also using proposition ^. 81 

ED H' [(•„.„ ., i = H S (3-"C) = (3-™) s W s (C) = 2-" = |C„, a „,..., ai r. We 

conclude that V is indeed a Vitali class for any subset E of C. 

Then, by Vitali's covering theorem (see PP theorem 1.10) for a fixed e > 0, 
there exists a finite or countable disjoint sequence of sets {Ui} from V such that 
either £ |C/ 4 | S = oo, or (J 17,) = 0, and also 

H s ( J B)<^|C/ l | s + £ . 

i 

Since the sets Ui are mutually disjoint and contained in C, and using (i) in 
lf2~TTJ) . it follows that 

mi* =J2n s (Ui) = n s (\J Ui) < n s (c) = 1. (2.12) 

i i i 

Therefore the other variant must be true: 
with U:= \J Ui 



H S (E\U) = 0. 



On the other hand 



H S (U\E) = H S (U) - H S (U n E) 

= H S {U)-{H S (E)-H S (E\U)) 

= n s {u)-n s {E)^^n s {u l )-n s {E) 

i 

i 

Also, observe that 

C,, ,, ., ••! " (c + f2a k 3 n - k ^j , (2.13) 

so by (I2~TUI) . x c „ an B1 = U~ n T l cp with Z = ££ =1 a fe 3"- fc . 

Therefore we see that all measurable sets E C K with H S (E) < oo are in 
the span of {U n T k ip \ n,k £ Z}. Since all integrable functions f E H can be 
approximated by simple functions, it follows that 

i7 = span {[/"TV | n, fc e Z} . (2.14) 
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3 Iterated Function Systems (IFS) and gap-filling 
wavelets 



The middle-third Cantor set C of Section 1 is a special case of an Iterated 
Function System (IFS). It falls in the subclass of the IFSs which are called 
affine. Specifically, let d G Z+, and let A be a d x d matrix of Z. Suppose 
that the eigenvalues \i of A satisfy |A»| > 1. Set N := |det A| . These matrices 
are called expansive. Then note that the quotient group Z d / A(Z d ) is of order 
N. A subset T> C Z d is said to represent the A-residues if the natural quotient 
mapping 

7: Z d -> Z d /A(Z d ) (3.1) 

restricts to a bijection 7 P of V onto 1 d /A(Jj d ). For example, if d = 1, and 
A = 3, then we may take either one of the two sets {0, 1, 2} or {0, 1, —1} as T>. 
The IFSs which we shall look at will be constructed from finite subsets S C Z d 
which represent the A-residues for some given expansive matrix A. If {A, S) is 
a pair with these properties, define the maps 

a s (x) := A _1 (x + s), s e S 7 x e M d . (3.2) 

Using a theorem of Hutchinson |15|. we conclude that there is a unique measure 
fi = A*(A,5) w hh compact support C = C(a,s) on ^ d such that 



se5 



or equivalently 



J f{x)dv{x) = ^yE / (3-4) 
The quotient mapping 

7: M d -> T d := R d /Z d (3.5) 

restricts to map C bijectively onto a compact subset of T d . The Hausdorff 
dimension h of /i and of the support C is 

k _ log #(5) 
log TV 

The system (C, /i) is called a Hutchinson pair, see lemma l3~51 

If d = 1, we will look at two examples: (i) (A, 5) = (3, {0, 2}) which is the 
middle-third Cantor set C in Section 1, and (ii) (A, S) = (4, {0, 2}) which is the 
corresponding construction, but starting with a subdivision of the unit interval 
/ into 4 parts, and in each step of the iteration omitting the second and the 
fourth quarter interval. As noted, then 

h (i) = lo S3( 2 ) = and tyii) = \> ( 3 ' 6 ) 
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for more details, see |TS] . 

Since the arguments from proposition 12.41 and 12.81 generalize, we will only 
sketch the general statements of results for the afhne IFSs, those based on pairs 
(A, S) in R d where the matrix A and the subset S C 7L d satisfy the stated 
conditions. The number h will be h = |°|j^ c ^jj ; i.e., the Hausdorff dimension 
of the measure /x, and its support C which are determined from the given pair 
(A, S). We will then be working with the corresponding Hausdorff measure 
H h , but now as a measure defined on subsets of M. d . The facts from Section 2 
apply also to this more general case in K. d , for example, property (1.2) for the 
middle-third Cantor set, now takes the following form 

AC = (J (C + s) (3.7) 

ses 

where AC := {Ax | x £ C}, and C + s := {x + s | x G C}, or equivalently 

C=\Ja s (C) (3.8) 

ses 

where a s (C) := {a s (x) | x £ C}. The conditions on the pair (A, S) guaran- 
tees that the sets in the union on the right-hand side in (|3.7|l or in (|3.8fl . are 
mutually non-overlapping. This amounts to the so-called open-set-condition of 
Hutchinson [Ej. The set TZ which is defined in Proposition l2.il in the special 
case of the middle-third Cantor set is now instead 

K=\J (J A~ n (C + k) = K= (J |J A~ n (C + k) (3.9) 

where C is the (unique) compact set determined by l|3.8|l . of Hutchinson's the- 
orem |15| . The properties of Proposition 12 . II carry over mutatis mutandis, for 
example, the argument from Section 2 shows that for every k £ 1 d and every 
n £ Z, 

K + A-' n k = n, and A n K = K. 

The Hilbert space H from Definition |2~SI is now H := L 2 (K,H h ). The unitary 
operators T and U from l2.8H2.9fl are now 

{T k f){x):=f(x-k), feH, x£K, k £ Z d (3.10) 

and 

(Uf)(x) = —L=f(A- x x), f£H,x£K. (3.11) 

v#( 6 ) 

The commutation relation from proposition 12 . 71 in its general form is 

UTklf- 1 = T Ak , k £ Z d . (3.12) 

We now need the familiar duality between the two groups Z d , and T d = M. d /Z d , 
which identifies points n £ Z d with monomials on T d as follows, 

z" = Z™ 1 Z% 2 ■ ■ ■ Z^ d — e i27rn i e i e i27rn 2 s 2 . . . e i"2itn d 6 d j-g -^2^ 
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Note that l|3.13[) identifies the torus T d with the d-cube 



{(9x,...,9 d ) |O<0;<1, i = l,...,d}. 

Since C is naturally identified with a subset of T d , we may view the monomials 
{z n | n £ Z d } as functions on C by restriction. We say that the system (A,S) 
is of orthogonal type if there is a subset T of 1 d such that the set of functions 
{z n | n £ T} is an orthonormal basis (ONB) in the Hilbcrt space L 2 (C, fi^ A s ^). 
If there is no subset T with this ONB-property we say that (A, S) is of non- 
orthogonal type. The authors of ^5] showed that (4, {0, 2}) is of orthogonal type, 
while (3, {0, 2}) is not. So for the Cantor set C 4 there is an ONB {z n \ n £ T} 
for a subset T of Z; in fact we may take 

{finite ^ 
Ke{0,i}j. (3.14) 

For the middle-third Cantor set C3 it can be checked that {z n | n £ Z} contains 
no more than two elements which are orthogonal in L 2 (C3,/x 3 ). 

Theorem 3.1. Let (^4,5) be an affine IFS in M. d , and suppose S has an exten- 
sion to a set of A-residues in 7h d . Let 

log #(5) 
1 log|det^| ; 

and let (C,/i) be as above; i.e., depending on (A,S), and let 1Z be defined from 
C in the usual way as in (|3.9|l . Assume further that 

Cn(C + fc) = 0, (k£Z d \{0}). (3.15) 

Then the system (A, S) is of orthogonal type if and only iff there is a subset T 
in Z d such that 

{(#(5))" /2 e l2 * A " k - x x c (A n x - t) I k £ T, (n = and I £ 7L d ) or 

(n> landers mod A for all s £ S)} (3.16) 
is an orthonormal basis in the Hilbert space L 2 (TZ,Ti. h ). 

Remark 3.2. The significance of the assumption l|3.15() is illustrated in \'26\. 
Also note that (|3.15|l is automatically satisfied if C = C(^4,6>) is contained in 
a 1 d -t\\e. This is the case for the example (^4,5) = (4, {0,2}), but there are 
examples in d = 2 where it is not. 

Proof. A simple check shows that 

TZ = |J {A- n {C + 1) I (n = and £ £ Z d ) 
or (n > 1 and £ ^ s mod A for all s £ S)} , 
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and the union is disjoint. Suppose {A,S) is of orthogonal type. We saw in 
Section 2 that the restriction of the Hausdorff measure Ti h to C agrees with the 
Hutchinson measure /i = fx^ A s \ on C = C^,s)- Hence density of {z n \ n e T} 
in L 2 (C, fi) implies density of {e^'-^ct^i I k £ T} in the subspace L 2 (C, U h ) 
of L 2 (lZ,H h ). Now the formula for 1Z implies that the functions in (|3.16|) are 

dense in L 2 (K,H h ). 

Suppose conversely that the family (|3.16|l is dense in L 2 (lZ,H h ). Then 
{z n | n € T} must be dense in L 2 (C, ji) since C is the support of Hutchinson's 
measure /i, and since fj, restricts 7i h . | 

Corollary 3.3. Let (C4,/i 4 ) be the Cantor construction in the unit interval 
I = T 1 defined by the IFS a (x) = f , o- 2 {x) = ^±2; i.e., by (A,S) = (4, {0,2}), 
and let 1Z be the subset of R defined in Ijff.ffl . Then the family of functions 

[2 n / 2 e i27r4nfa x c (4"a; - I) \ k € {0, 1, 4, 5, 16, 17, • • • }, 

£ f Z if n = 1 , . 

fe \ Z\ (4Z + {0,2» £/ »>1. J 1 > 

forms an orthonormal basis in the Hilbert space 

Proof. This is a direct application of the theorem as the subset 

T = {0,1,4,5, 16, 17, •••} 

from (|3.14|l and (|3.17|) satisfies the basis property for C4, ^ 4 by Theorem 3.4 in 

Mm 

The next result makes clear the notion of gap- filling wavelets in the context of 
iterated function systems (IFS). While it is stated just for a particular example, 
the idea carries over to general IFSs. Note that in the system (|3.18ll below of 
wavelet functions, the two ip 2 an( i ^3 are gap-filling. 

Corollary 3.4. Let C = C4 be the Cantor set determined from the IFS, o~o{x) = 
j, o~2(%) = ^4^> from the previous corollary. Then the three functions 

M*) = =Xc(4z)-Xc(4^-2) (3.18) 

iP 3 {x) : =V2 Xc {4x-3) 

generate an orthonormal wavelet basis in the Hilbert space L 2 (Ti,m). Specifi- 
cally, the family 

{2* V 4 (4 fe a; - I) I i = 1, 2, 3, k € Z, I € z} (3.19) 
is an orthonormal basis in L 2 (TZ,Ti.^). 
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Proof. We noted that our results in propositions 12.41 and 12.81 apply more gen- 
erally to IFSs of affine type. So the result amounts to checking the general or- 
thogonality relations for the functions mo, mi, m,2, TO3 on T which define wavelet 
filters for the system in (|3.19|) . Note that from (|3.19fl the subband filters {mi}f_ 
are as follows, zeT: 



m (z) = 

toi(z) = 

m 2 (z) = 
to 3 (z) = 

Since the 4x4 matrix in the system 



^ <1 + "» 



/ m (z) \ 

TOl(z) 
TO 2 (z) 

\ m 3 (2) / 





V o 



z . 



1 


M 






1 \ 











z 












z 2 





i y 






z 3 J 



is clearly unitary, the result follows from a direct computation; see also the proof 
of theorem 16. 21 

To verify that the Ruelle operator R — R mo given by 



- \m (w)\ 2 f(w) 



E 



w 2 + w 2 



satisfies the two conditions 

(a) dim{/ S C(T) \Rf = f} = l, and 

(b) for all A e C, |A| = 1, and A ^ 1, dim{/ e C(T) | Rf = A/} = 0, we 
may again apply the theorem from |14| or the results of section 6 below. | 

For the more general affine IFSs the results above extend as follows. 
Consider the affine IFS {o~i) p i=1 with 



(x e 



where N > 2 is an integer and (ai)f=i are distinct integers in {0, ...,N — 1}. 
Then by , there is a unique compact subset if of K which is the attractor of 
the IFS, i.e., 

C = Uf =1 a,(C). 

Actually, one can give a more explicit description of this attractor, namely 
C = {J2djN- j | dj e {ai, a p }, j > 1}. 
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Since the digits aj are distinct and less then TV, K is contained in [0, 1], and the 
sets Oi{K) are almost disjoint (they have at most one point in common, those 
of the form k/N for some k G {1, TV — 1}. 

The Hausdorff dimension of K is log^ p. 

Now consider the set 

1Z = { djN^i | to e l,,dj G {ai, a p } for all but finitely many indices j} 

j>-m 

1Z is invariant under integer translations 

!Z + k = lZ, (k e Z), 
and it is invariant under dilation by N 

Nil = 11. 

Endow 1Z with the Hausdorff measure U s for a = \og N p, and on L 2 (K,H S ), 
define the translation operator 

Tf(x) = f{x - 1), (x ellje L 2 (1Z,H S )), 
and the dilation operator 

uf(x) = ^f (£) , (x e nj g l 2 (iz,h s ))- 

These are unitary operators satisfying the commutation relation 

UTU^ 1 = T N . 

Let </? := Xc- The function p is an orthogonal scaling function for L 2 (1Z, H s ), 
with filter 



mo(z) = J-yz a \ (3.20) 
VP'S 

so it satisfies the following conditions: 

1. [Orthogonality] 

(T k (p | if) = S k , (k € Z). 

2. [Scaling equation] 



3. [Cyclicity] 

span{U~ n T k ip \n,keZ} = L 2 (K, H s ). 
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Next, we define the wavelets. For this, we need the "high-pass" filters 
mi, ...,tojv-i such that the matrix 

* (rn^))^, 



is unitary for almost every z. (p = e 2m / N ). 

First, we define the filters for the gap-filling wavelets -01; tpN-p- The set 
G = {0, N — 1} \ {cti, a p } has N — p elements. We label the functions 
z i— > z d for d £ G, by mi, rriN-p- 

The remaining p— 1 filters are for the detail-filling wavelets. Let rj = e 27 ™/ p . 
Define 



^E^- 1 ^* G {1, — ,P- 1})- 



We have to check that 

l^m.W^H^,, (*GT,i,jG {0,...,N-1}). (3.21) 

1D JV =Z 

For this we use the following identity: 

^2 w k = 0, (zeT,^0 mod AT). 



Therefore, if fi(z) = h = Eilo 1 ft 2 *' then 

JV-l iV-l 

Applying these to the filters m*, (i € {0, A" — 1}), we obtain (|3.21() . 
With these filters, we construct the wavelets in the usual way: 

i>% = U~ 1 m i (T)ip, (ie{l,...,AT-l}), 

and 

{U m T n ^ l | m,n G Z, i G {1, .... JV - 1}} 

is an orthonormal basis for L 2 (TZ,H S ). 

Let AT G Z + be as above, and consider S = {ai, • • • , a p } C {0, 1, 2, • • • , AT — 
1}. A second subset B = {b\, ■ ■ ■ , b p } c Z is an A^-dual if the pxp matrix 

M^) = ^(cx P (^^ (3.22) 



1<j./c<P 



is unitary. When A" and 5 are given as specified, it is not always true that there 
is a subset B c Z for which M^(S,B) is unitary. If for example AT = 3 and 
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S = {0,2}, then no B exists, while for N = 4 and S — {0,2}, we may take 
B = {0, 1}, and 

is of course unitary. 

Lemma 3.5. \16lj Let N and S be as specified above, and suppose 

B = {&!,•■■ ,6 P } C Z 
is an N-dual subset. Suppose £ B, and set 

{finite ^ 
mN* \rneBj. (3.23) 

Let (C,/i) = (C(jv,5), M(at 5)) ^ e ^ e Hutchinson pair. Then the set of functions 
{z n \ n £ A} is orthogonal in L 2 (C,/i); i.e., 

« n -"'d/x(*) = * n y 1 n,n'eA (3.24) 

where we identify C as a subset of T 1 via 

C 9 6 — ► e i27re € T 1 . 
Proof. Set e(0) = e l27Tf> , and for fc € E 

B(jfe) := J e(k0)dn{6). (3.25) 
c 

Using l]3.4[l. we get 



^=^°UJ fl UJ' (3 - 26) 

where mo is defined in l|3.20fl . 

If n, n' € A, and n ^ n', we get the representation 

n' - n = b' - b + mN e , b, b' e B, m, I e Z, i > 1. 
As a result, the inner product in L 2 (C,/u) is 

1 (V -b\ fri-n\ 

Since the matrix Mn(S,B) is unitary, 



when &' ^ 6 in B, and the result follows. 
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Even if the matrix Mn(S,B) is unitary, the orthogonal functions {z n | n G 
A} might not form a basis for L 2 (C,[i). From |16|. we know that it is an 
orthonormal basis (ONB) if and only if 

\B(£-n)\ 2 = 1 a.e. £ G R. (3.28) 

nGA 

Introducing the function 



* beB 



m--= l -Y,\ B ^- n )\ 2 ' ( 3 - 29 ) 

and the dual Ruelle operator 



P beB 



m 



N 



g-6 
N 



(3.30) 



we easily verify that f2 and the constant function 1 both solve the eigenvalue 
problem Re(f) = f, both functions f2 and 1 are continuous on R, even analytic. 

Theorem 3.6. If the space 

{/ G Lip(R) I / > 0, /(0) = 1, R B (f)=f} (3.31) 

is one-dimensional, then A(= Ajy(B)) induces an ONB; i.e., {z n \ n G A} is an 
OAS in L 2 (C,/i). 

Proof. The result follows from the discussion and the added observation that 
O(0) = 1. This normalization holds since £ B was a,ssiim.6d, and so (co | Cn) w — 
for all n G A\{0}. | 

Definition 3.7. A S-cycle is a finite set {zi, Z2, ■ ■ ■ , Zk+i} C T, with a pairing 
of points in B, say bi, &2j • • • , 0/c+i € such that 

Zi = o-- bi (z i+ i), z k +i=zi, (3.32) 

and |mo(2 : i)| 2 = p. Equivalently, a B-cycle may be given by . . . , C R 
satisfying 

ee bi + N^ modiVZ 
(N k - 1) = + + • • • + AT* -1 **! mod AT fe Z. 

Theorem 3.8. Let N G Z + , N > 2 be given. Let S C {0, 1, • • • , N - 1}, and 

suppose there is a B C Z swc/i i/iai € i3, #(<S) = #(<6) = P, and the matrix 

2nab\ 



is unitary. Then {z n \ n G Ajv(i3)} is an CWZ? /or L 2 (C,/i) where An(B) is 
defined in i/ i/ie onZj/ B-cycles are the singleton {1} C T. 
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Proof. By Theorem l3.6l we need only verify that the absence of 23-cycles of order 
> 2 implies that the Perron-Frobenius eigenspace (|3.30|l is one-dimensional. But 
this follows from 6 Theorem 5.5.4]. In fact, the argument from Chapter 5 in 
|H] shows that the absence of K-cycles of order > 2 implies that the 2?-Ruelle 
operator R$ with cr_h(^) := 

(ife/XO = -E (\mo(o-- b (0)\ 2 f(v-b(0) 

P beB 

satisfies the two Perron-Frobenius properties: 

(i) the only bounded continuous solutions / to Rs(f) = f are the multiples of 
1, and 

(ii) for all A 6 T\{1}, the eigenvalue problem Rs(f) — A/ has no non-zero 
bounded continuous solutions. | 

Example 3.9. (An Application) Let N — A, S = {0, 2}, and B = {0, 1}. Then 
A 4 (B) = {0,1,4, 5, 16, 17, 20,21, •••}, and 



{Rem) = cos^(2<) / I J ] W(2<) / 

and there is only on S-cycle, the singleton {1} C T. Recall from ^5] that 
the Hutchinson construction of (C, ^) identifies C as the Cantor set arising by 
the subdividing algorithm starting with the unit interval / dividing into four 
equal subintervals and dropping the second and the fourth at each step in the 
algorithm. The measure n is the restriction of H.? to C, and it follows from 
the last theorem that {z n \ n G A 4 (S)} is an ONB for L 2 (C,^). The dual 
system {cr_{, | b £ B}; i.e., er (£) = f j — generates a Cantor subset 

Cg C [—1,0] also of Hausdorff dimension i. Note that the fractional version of 
the Ruelle operator i?g does not map 1-periodic functions into themselves; in 
general 

(ife/XO? 6 (%/)(£ + !); 



in fact 



so 



(RbM + 1) = cos 2 (2<)/ (j^J + sm 2 (2<)/ 



holds only if 

cos 4 (2<)/(^) =sin 4 (2^)/ 



4/n ,w/£ + M „.-„4,o_mWC-1 



The following tables of similar examples is included hopefully offering the 
reader a glimpse of the variety of examples, all of orthogonal type. The tables 
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also offers some insight into the duality between the two systems, one in the 
a;- variable and the other in the Fourier dual variable £ 

N p S B M N (S,B) Hausdorff Dim. 

1 1 



i I x x i I 

2 




4 2 {0,2} {0,1} 

6 2 {0,3} {0,1} ^jlj 1 1 I log 6 (2) 

6 2 {0,1} {0,3} -±=( J l ,) log 6 (2) 



6 3 {0,2,4} {0,1,2} ^ « i log6 ( 3) 

\ 1 S3 S3 / 

where ( 3 = exp 

N p A N (B) 

4 2 {0,1, 4,5, 16, 17,20, 21, •••} 

6 2 {0,1,6,7,36,37,42,43,- • •} 

6 2 {0,3,6,9,36,39,42,45,- • ■} 

6 3 {0, 1,2, 6, 7, 8, 36, 37, 38, 42, 43,44, •••} 

<)/(!) +sin 2 (27r£)/(^i) 
^)/(l)+-n 2 (^)/(^) 
f)/(|) + sin 2 (f)/(i^) 

^(¥)/(¥) + ^(¥)/(¥) 

where W(0 := 2cos2(2 ^ ) 3 " sin2(3 ^ ) . 



iV 


P 




4 


2 


cos 2 


6 


2 


cos' 


6 


2 


cos : 


6 


3 





4 A Generalized Zak- Transform 

The notion of filter is imported into math from signal processing. It has now 
been well adapted to wavelet analysis: In the familiar dyadic case the two 
wavelet functions ip (the father function) and ip (the mother function) are 
in the Hilbert space L 2 (R). In the TV-adic case, the wavelet functions are 
tp, ipi, ■ ■ ■ , tpN-ii an d there are known conditions for when these functions are 
in L 2 (M). The starting point is the scaling identity 1)1. 6fl satisfied by ip. Intro- 
ducing the wavelet filter m${z) — ^^akz k as a function on T = R/27rZ, and 

the transfer operator R mo in l|1.10|l . we note that necessary conditions for (p to 
be in L 2 (R) are |rao(l)| = %/iV, the low-pass condition, and i? mo (i) = 1, where 
1 is the constant function 1. Let 5± denote the Dirac mass at z = 1. It follows 
that 5iR mo = S\. But what if, for some mo, Si does not satisfy this, so called 
low-pass condition; but rather there is some other probability measure v on T 
which is R mo -invariant; i.e., satisfies vR ma = v i an d which is singular with full 
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support. A main point in our paper is that this alternative introduces fractal 
analysis into the wavelet construction. 

Traditionally, the Zak-transform is a standard tool of analysis in L 2 (K), 
and in this section it is extended to the abstract case of Hilbert space. If 
{T k : k G 1} denotes translation (T k f)(x) = f(x - k), x G R, k G Z, / G L 2 (R), 
we set 

(Zf){z,x) = Y,* k Tkf(x), z&T,xGl, 
fcez 

and we check that Z defines a unitary isomorphism of L 2 (R) onto L 2 (T x 7) 
where T is the torus, and I the unit-interval I = [0, 1). The measure on T is 
Haar measure, denoted fi; i.e., 

J ■ ■ ■ dfi(z) = J ■ ■ ■ dO, where z = e l6 . 

T 

Let Ti. be a Hilbert space, U a unitary operator in Tt, and T: Z — > U(TL) a 
unitary representation. Let N >2, and suppose that 

(77W- 1 = TV fc , fceZ. (4.1) 

In general, the inner product in a Hilbert space 7i will be written (■ | ■). If 
f G Ti, then the form g — * (/ | g) is taken linear on 7i; and (/ | /) = ||/|| . 
For /j € H, i = 1, 2, we introduce the following function 

p(A,/ 2 )(z) |/ 2 ) (4.2) 

defined formally for z G T. Let mo G L°°(T) be given, and suppose that 

i^|m (w;)| 2 = l, a.e. zeT. (4.3) 

UJ N =Z 

Note that the sum in l|4.3|l is finite, since for each z G T, the equation w N = z 
has precisely solutions. In fact, the cyclic group Zn acts transitively on this 
set of solutions {w}. 

(We will work with the torus T in anyone of its three familiar incarnations: 
(i) {z £ C | \z\ — 1}, (ii) the quotient group R/27rZ, or (iii) the period interval 
[0, 27r) via the identification z = e l6 ' . With this identification, we get the familiar 
description of the set {w G T | w N — z(— e l9 )} in the summation on the left- 
hand side of as the N distinct frequency bands { e+ ^ rk |fc = 0,l,---,A r -l} 
The sub-interval [0, ^) represents the low frequency band.) 

The operator mo(T) is defined from the spectral theorem in the usual way: 
If the spectral measure of T is denoted Ej>, then Ex is a projection valued 
measure on T, and we have the following three identities: 

||/|| 2 = |||£; T (dz)/|| 2 , fen, (4.4) 

T 
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T k = J z k E T (dz), keZ, (4.5) 

T 

and by functional calculus, 

m (T) = J m (z)E T (dz) (4.6) 

T 

If _ 

m Q (z) = J2a k z k (4.7) 

fcez 

is the Fourier series of too, it follows that toq(T) = J2k£Z a kT k is then well 
defined. 

The following operator i? = i? mo , called the Ruelle operator, (see 1S 
acting on functions h or T as follows, 

(ffl)(z) = i^KW| 2 %), ^T. (4.8) 

Let 1 denote the constant function 1 on T. Then condition (|4.3[1 amounts to 
the eigenvalue equation 

R(l) = i (4.9) 
On the Hilbert space TL, we introduce the operator 

M := U^moiT). (4.10) 

It is called the cascade approximation operator. In the special case when T k f(x) — 
f(x - k), (17/) (x) = i\T */(#), and H = L 2 (R), then 

(Mf)(x) = VNj2a k f(Nx-k) (4.11) 

k 

where {a k : k € Z} is the sequence of Fourier coefficients of mo, see (|4.7(l . In this 
case M is also called the wavelet subdivision operator. The following general 
lemma applies also to the case of wavelets in L 2 (1R). The advantage of l|4.10|l 
over H4.11fl is that l|4.10|l is defined for all systems U, T satisfying (|4.1|) and 
applies in particular to our present fractal examples. 

For measurable functions £ and rj on T, formula (|4.f>|) represents the usual 
functional calculus; i.e., £(T) = f T £(z)ET(dz). Setting 7r(£) := £(T), we get 
^i^v) = n (0 7T ( r l)^ tt(1) = I = the identity operator, 7r(£) = 7r(£)* = the 
adjoint operator. These properties together state that 7r(= tte t ) defines a 
^-representation of _L°°(T) acting on the Hilbert space Ht of the translation 
operators {T k : k G Z}. 
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Lemma 4.1. Let TL, T , U , and mo be as described above, and let the operators 
M and R be the corresponding operators; i.e., the cascade operator, and Ruelle 
operator, respectively. Then the identity 

R(p(fuf2))=p(Mf u Mf 2 ) (4.12) 

holds for all /1./2 € Ti., where the two sides in are viewed as functions 

on T. 

Proof. Let £ G C(T). Then it follows from JO} and JOJ that 

J &)p(h,h){zW{z) = (fi I e(T)/ 2 ) (4.13) 

T 

where /i is the Haar measure on T, (• | •) is the inner product of Ti, and £(T) = 
J T ^(z)dET(z). Using this, in combination with (|4.1U|) . we therefore get 

p{Mf u Mh){z) = ^z k {T k U^ 1 m (T)fi | U~ 1 m (T)f 2 ) 
= Y, zk i T Nkm n {T)h I m (T)/ 2 ) 

feGZ 

= ^ E p(A-KI 2 (t)/ 2 )h 

w GT 

jv 

w = z 

and therefore 

/" Z(z)p(Mh,Mf 2 ){z)dn{z) 

T 

= J (i(z N )p(fu\™o\ 2 (T)f2)(zW(z) 

T 

= /" ^ImoOaOlV/i./aX^MCz) 

T 

= Ji(z)R(p(fiJ 2 ))(z)dn(z). 

T 

Since this is valid for all £ G C(T), a comparison of the two sides in the last 
formula, now yields the desired identity l|4.13l) . | 

Remark 4.2. An immediate consequence of the lemma is that if some ip G Ti 
satisfies Mip = ip, or equivalently 

u^ = Y. akTk v ( 4 - 14 ) 

feez 
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then the corresponding function h :— p(<p, <p) on T satisfies R(h) — h. Recall 
that / 14.14J) is the scaling equation. If further the functions {T^ip : k € Z} on 
the right hand side in (|4.14|l can be chosen orthogonal, then 

p{v>,<f){z) = {<? I <p)n = WfWn ( 4 - 15 ) 

is the constant function, and so we are back to the special normalization condi- 
tion 14.3)) above. 

The function p((p, ip) is called the auto- correlation function since its Fourier 
coefficients 

J z- k p(<p, tp)(z)dii(z) = (T k <p | ip) H (4.16) 

T 

are the auto-correlation numbers. 

Lemma 4.3. Let mo S C(T) be given, and suppose holds, (a) Then there 

is a probability measure v = v mo depending on mo such that 

J \dv = J R(£)dv (4.17) 

T T 

holds for all£e C(T). 

(b) If mo is further assumed to be in the Lipschitz space Lip(T), then the 
following limit exists 

-. n— 1 

lim - V fj,R k = v (4.18) 

n — >oo 77, ' J 

where ,R^) := = /* » ^ Foot ~e, and <fte c n- 

— T 

vergence in is in the Hausdorff metric (details below). The measure v 

satisfies J^.J7| /. Moreover, when then operator R from C (T) to C (T) has Perron- 
Frobenius spectrum (i.e., 1 is the only eigenvalue of absolute value 1) then the 
limit 

lim (j,R n = v, (4.19) 

n — >oo 

exists and gives the unique invariant measure v . 
Definition 4.4. A function £ on T is said to be in Lip(T) if 

£(e is ) - £(e rf ) 



Dt; := sup 

— 7T<S<i<7T 



s - i 

The Lipschitz-norm is ||£|| Li := D£ + $(1). 



< oo. (4.20) 
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The Hausdorff distance between two real valued measures Vx,v% on T is 
denned as 

distjj aus (^i, v%) — sup < J ^dv\ — J ^dv2 | £ G Lip(T), £ real valued, and D£ < 1 
It' 

(4.21) 

Hence the conclusion is the lemma states that if too £ Lip satisfies (|4.3J) . then 

lim distTT ail „ (>,^iT) = 0. (4.22) 

Note that both the Ruelle operator R and the measure v depend on the func- 
tion Too- Condition (|4.3|) states that the constant function 1 is a right-Perron 
Frobenius eigenvector, and l|4.17|l that v is a left-Perron-Frobenius eigenvector. 



Proof. (Lemma The lemma is essentially a special case of the Perron- 
Frobenius- Ruelle theorem, see Also note that an immediate consequence 
of H4.17fl is the invariance 

' \{z N )dv{z)= h(z)dv{z) (4.23) 



A key step in the proof is the following estimate: Let ||£|| :— sup z6T |£(z)|; 
then the following estimate 

||^ll Lip <^rlieil Lip + 2i?(|TOo| 2 )-|ieil (4.24) 

hold for all £ £ Lip(T) and all n G Z+. We leave the details of the verification 
of l|4.24|l to the reader; see also [Hj, and [7|. | 

For functions to and m! on T, define the form (to, to') as a function on T as 
follows 

(to,to')(z) :— — m(w)m'(w). (4-25) 

idGT 

W N =z 

Lemma 4.5. Let N > 2, and let (Ti, U, T) be a system which satisfies the 
commutation relation ■ Let to and to' be functions on T which both satisfy 
the normalization condition Let M — M m be the cascade approximation 

operator, M = U^ 1 m(T), and suppose the scaling identity Mip = ip has a non- 
zero solution in H such that the vectors Tj~tp, k G Z. are mutually orthogonal. 
Let M' = M m > = U- l m'(T). 
Then 

p( ( p,M' ( p)(z) = \\<pf n -(m,m!)(z), z e T. (4.26) 
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Proof. Using Mip = ip, and p(<p, <p>){z) = ||y||^, we get 
p(<p,M'<p)(z) = p(M(p,M'(p) 

= "ft 51 Tn[w)m'(w)p{ip,(f)(w) 

W N =Z 



= WvWn'Tf H m(w)m'(w) = \\cp^- (m,m')(z) 



which is the desired identity (|4.26fl in the conclusion of the lemma. | 
Corollary 4.6. With the assumptions in Lemma \4-5[ set 

m { - n \z) := m{z)m{z N ) . . . m^"^) (4.27) 

and 

m l{n \z) := m'{z)m\z N ) . . . m'^^). 

Then 

p(<p,M' n <p)(z) = M 2 H ~ ^Hm' {n \w) (4.28) 

w Nn =z 

Proof. A direct iteration of the argument of Lemma 14.51 immediately yields the 
desired identity (|4.28j) . | 

Example 4.7. Let N = 3, and let the functions m and mf be given by m(z) = 
^| and m'(z) = z 3 m(z). Let 

H = L 2 {1l,(dx) s ), 3 = log 3 (2), (4.29) 

and ip = Xc! i- e -j V € W is the indicator function of the middle-third Cantor set, 
Cel. Let Tf{x) = f(x - 1), and (Uf)(x) = ^/(f )■ Then the conditions in 
Lemma ETSI and Corollarv k. Bl are satisfied for this system. Specifically \\<p\\ n = 1, 
Mip = (p where M = [/ _1 m(T); i.e., 

<p(x) = ip(3x)+ip(3x-2), (4.30) 

or 

3C = CU(C + 2). (4.31) 

Also notice that 

M' = U- x m'{T) = TM. (4.32) 

It follows that 

p(ip,M' n ip)(z) = z W+^ 2 +-+^- 1 = z^r 1 - 

and therefore 

(ip,M' n ip) = / p(ip,M' n ip)(z)dn(z) = 
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for all n 6 Z + . In contrast to the standard cascade approximation for L 2 (R) 
with Lebesgue measure, we see that the cascade iteration on the Cantor function 

Xci Le -> 

Xc ,M' Xc ,M' 2 Xc ,... (4.33) 

does not converge in the Hilbcrt space Ti s = L 2 (1Z, (dx) s ). In fact the vectors 
in the sequence (|4.33|l are mutually orthogonal. 

While the measure v € M(T) from Lemma [4.51 is generally not absolutely 
continuous with respect to the Haar measure /j, on T, the next result shows that 
it is the limit of the measures |to^"^(z)| dfi(z) as n — > oo if the function m is 
given to satisfy (|4.3|l and if the sequence is defined by (|4.27|l . 

Proposition 4.8. Let m e L°°(T) be given. Suppose \4-3[ holds, and let v be 
the Perron- Frobenius measure of Lemma \4-.5\ i.e., the measure v — v m arising 
as a limit \4.19^ . Then 

lim / £(z) m {n) {z) 2 dfi[z) = [ t{z)dv{z) (4.34) 



holds for all£e C(T). 

Proof. Calculating the integrals on the left-hand side in (|4.34|) . we get 
J £ dfJ,= J £iT"(i)d/i = J R n {£)dli—^ J idv 

T T T T 

where (|4.18|) was used in the last step. This is the desired conclusion (14.34ft of 
the proposition. | 

Corollary 4.9. If m(z) = ~|r * s ^ e function from Example \4- ?1 then the 
substitution z — e lt yields the following limit formula for the corresponding 
Perron- Frobenius measure v , written in multiplicative notation: 

1 " 

lim 7T- TT (1 + cos ( 2 ' 3 fc *)) = dv(t) (4.35) 

fe=l 

Proof. The expression on the left-hand side in (|4.35(l is called a Riesz-product, 
and it belongs to a wider family of examples; see for example [7], and [5]. 
The limit measure v is known to be singular. It follows, for example from 
jS]. A computation shows that the Fourier coefficients v(n) := L z n dv{z), 
are real valued, satisfy v(0) = 1, v(X) = 0, £>(— n) = v(n), P(3n) — v(n), 
v(Zn — 2) = ^(n), z>(3n + 2) = \v{n), for all neZ. Hence 

1 - 

lim |i>(fc)| 2 = = sum of atoms; (4.36) 

k— — n 
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i .e., W({z})\ 2 = 0. The last conclusion is from Wiener's theorem, and implies 
that v has no atoms; i.e., v({z}) = for all z E T. 

To check <|4.36ll some computations are required. Let 

s n =|£>(0)| 2 + |/>(l)| 2 + ... + |*>(n)| 2 , (neN). 

Then, using the recursive relations for we have 

3 „ + i 

s 3 „ +1 = S3 „ + \Hk)\ 2 

fc=3 n + l,fc=0 mod 3 

311+1 gn+l 

+ E i^ fc )i 2 + E i^)! 2 

fe=3"+l,fe=2 mod 3 fe=3"+l,fe=-2 mod 3 

115 

< S 3 „ + S 3 n + -S 3 n + -S 3 n = -S 3 „. 

By induction 

■»-<(§) »o. 

Now take fc arbitrary then for some n, fc is inbetween 3 n and 3™ +1 so 

Sk < £3^ < (|)" +1 ao = /5\" 5^ 
jfe _ 3" _ 3™ 1 6/ 2 °' 



which shows that Sk/k converges to and this proves 

In summary, the measure v is singular and non-atomic. In the next section 
we show that v has full support. 

Moreover there are supporting sets for v which have zero Haar measure as 
subsets of T. Concrete constructions are given below. 

It follows from the recursive relations for the numbers v{n) that P(2k-\-X) = 0, 
k G Z; i.e., that all the odd Fourier coefficients vanish. 

Each integer n G Z + has a representation of the following form: 

n = l + h ■ 3 + h ■ 3 2 + • • • + l p ■ 3 P , U E {0, -2, 2}, i<p,l p e {1, 2}. (4.37) 

The representation is not unique (because for example 1 = 3-1 — 2) but unique- 
ness is obtain if we impose that l v -\ ^ — 2 when l p = 1. We define a counting 
function #(n) which records the occurrence of values —2 and 2 for the 'trigets' 
li. Hence, if n G Z + is even, then 

v{n) = 2-* {n) . 

Now, introduce the following sequence of functions 

5fe (z):=z 2 - 3fc -i, zel (4.38) 
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with inner products as follows with respect to the measure v: 



(9k | 9i) v = I g k (z)gi{z)dv{z) = -5 k ,i (4.39) 



T 

It follows from this that the series 



^) : =EfcX7^W, zeT, (4.40) 

fe=0 

is convergent in L 2 (T, v) with 



°° t 2 

||2 1 n n 2 7T 



Using the Riesz-Fisher theorem, we get a Borel subset A cT, v(A) — 1; i.e., 
^(T\v4) = 0, and a subsequence n\ < n,2 < < . . . , rii — > oo, such that the 
series 

"> i 

fe=0 

is pointwise convergent, i — * oo, for all z e A But note that 

y J_ 9fe (z) = y -J- z 2 -3 fc _ ly _JL . (4 . 4 3) 

Using now Carleson's theorem about Fourier series on T with respect to Haar 
measure fi (=Lebesgue measure), we conclude that there is a Borel subset B C 
T, fi(B) = 1; i.e., n(T\B) = 0, such that the series J2k=o Fpj" z2 3 i s pointwise 
convergent, n — > oo, for all z £ B. But identity (|4.43() implies that A n B = 
0, and so fJ,(A) = 0. The supporting set A for the measure v has Lebesgue 
measure zero; and moreover the two measures v and \x (=Haar measure on T) 
are mutually singular. | 

It can be shown, using a theorem of Nussbaum that the Ruelle operator 

has Perron- Frobenius spectrum on C(T), specifically that 

E R (1) = {fe C(T) | Rf = /} = Ci, 
and if A G C, |A| = 1, and A ^ 1, that 

Sa(A) = {/ S C(T) | Rf = A/} = 0. 
As a consequence we get that 

ll^ l (/)-K/)ll Lip <^ll/ll Lip 

holds for all / 6 Lip(T) where ||-|| Li p denotes the Lipschitz norm on functions 
on T. 
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5 The support of Perron- Frobenius measures 

In this section we consider the support of the measure v. Caution: By the 
support of v, we mean the support of v when it is viewed as a distribution; 
i.e., the support of v is the complement of the union of the open subsets in T 
where v acts as the zero distribution, or the zero Radon measure. Even though 
the support of v may be all of T, there can still be Borel subsets E C T with 
v(E) = 1, but E having zero Lebesgue measure. 

Theorem 5.1. Let mo € Lip(T). Suppose \4-&t holds and mo has finitely many 
zeros and let v be the Perron- Frobenius measure of Lemma \4-'A Then exactly 
one of the following affirmations is true: 

(i) The support of v is T. 

(ii) v is atomic and the support of is is a union of cycles C = {z%, . . . , z p } with 
Zi = z 2 , ■ ■ ■ , Zp_ x = z p , z^ = z\, and \m \ 2 (zt) = N for all i £ {l,...,p}. 
(Such cycles are called (mo, N)- cycles). 

Proof. To simplify the notation, we use 

W:=\m Q \ 2 . 

oo 

Note that proposition 14.81 gives us v as an infinite product j^J W(z N )dfi. 

fc=i 

In the next lemma we analyze the measures given by the tails of this product. 
Lemma 5.2. Fix n > 0. We then have: 

(i) For all f G C(T) the following limit exists, and defines a measure on T: 
v n {f) := hm^oo / T W(z N ") ■ ■ ■ W(z Nn+k )f(z)dv 

(ii) v n (f) = v(R™f),(f G C(T) where R\ is given as in Q4.8H but with mo = 1. 

(iii) J T f(z)W(z) ■ ■ ■W(z N ^ 1 )dis n = f T f(z)dv, f £ C(T). (So v is absolutely 
continuous with respect to v n .) 

(iv) IinWooM/)=Ai(/),/GC(T). 

Proof of Lemma. We can use a change of variable to compute 
hm J W (z r ) ■ • • W f(z)dfi 

= lim / W(z)---w(z N ")R^f(z)dtJ, 
= v(R»f). 

This proves (i) and (ii). (iii) is immediate from proposition ^. 81 
For (iv) note that for / e C(T) and neff. 

^)^E'z(^) 

fe=0 v 7 

therefore Rff converges uniformly to L fdfi. 
Then, with (ii), lim B -.ooi'fl(/) = I 
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We continue now the proof of the theorem. We distinguish two cases: 
Case I: All measures v n are absolutely continuous with respect to v. In this 
case we prove that the support of v is T. Assume the contrary. Then there is 
an open set U with v(U) = 0. This implies v n (U) = for all n. 

Take / S C(T) with support contained in U. Then v n (f) = 0. Take the limit 
and use lemma I^^T ivl . it follows that /«(/) =0. As / is arbitrary, /j,(U) = 0. 
But this implies U — 0, so the support of v is indeed T. 

Case II: There is an n € N such that v n is not absolutely continuous with 
respect to v. This means that there is a Borel set E with v{E) = and 
v n {E) > 0. 

We prove that there is a zero of W^ n \ call it Zq, such that v n ({za}) > 0. 
Suppose not. Then, take E' = £\zeros(WW), v{E') = u n {E') > 0. 

The measure v is regular so there is a compact subset K of E' such that 
v n {K) > 0. Of course v(K) = 0. Since K has no zeros of W^ n ' and this is 
continuous, is bounded away from on K. Then, with lemma l5"2T iii) 

= / wm* = I w{n){z) ^h) dVn = Vn{K) 

K K 

which is a contradiction. 

Thus, there is a z S zeros(W( n )) with v({z }) > 0. 

We know also that v(f) = v(Rf) for all / S C(T). By approximation 
(Lusin's theorem) the same equality is true for all bounded Borel functions. 
Then 

o < v(x {zo} ) = v{Rx {za} ) = v ( ^ w Hx {zo }H ) = ] ^j^v(x {z »})- 

\ W N =Z ) 

(5.1) 

Therefore W(z ) > and v({z$}) > 0. By induction, W(ztf") > and 
v{{zf}) for all k e N. 

Since (|4.3|) holds, W(z) < N for all z € T; so from the previous computation 
we obtain 

Since z/ is a finite measure, the orbit ^Zq H \ k € n| has to be finite so the 
points zq, Zq , . . . , Zq V will form a cycle, Zq P+ = zq. Also, 

„({«>}) = - ({^ P+1 }) > - }) >••■>- ({*}) 

hence all inequalitites are in fact equalities and with Q5.ip . this shows that 

W'(*^*)=.2V'far*e{0,...,p}. 

We are now in the "classical" case and we can use corollary 2.18 in [H] (See 
also to conclude that v must be atomic and supported on cycles as mentioned 
in the theorem. | 
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Lemma 5.3. Let mo,m' be Lipschitz, with finitely many zeros. Suppose 



-Rmo(l) — 1 — Rm' (l), 

and suppose there are no mo or mQ-cycles. Assume in addition that mo and m' 
have the same Perron- Frobenius measure v, then \m$\ = \m' \. 

Proof. With W := |too| 2 , we have, from proposition ^. 81 for / G C(T): 
lim / f{z)W {n) {z)d^{z) = { f{z)dv{z). 



Then 



Jim / f(z)W {n) (z N )dfi = Jirn / R 1 f{z)W {n) {z)d^ 

T 

Rxf{z)dv. 



Take / G C(T) such that / is zero in a neighborhood of zeros(mo). Then ^ ^ 
is continuous. So 

lim 



l^W^(z)d f i(z)= f f{z) 2 d V (z) 



\m (z)\ 



|m (z)| z 



On the other hand 



lim 

n — >oo 



rr^^WW = lim 

|m (z)| ™^°< 



/(*) 



|ra (z)| z • • • mofz* ) 



m (z)\ 

= lim / f{z)W {n - 1) {z N )dv{z) 



Rif(z)dv. 



dn 



Thus 



/(*) ^ 
|m (z)|" 



dviz) = / R 1 f{z)dv{z) 



(5.2) 



for all / G C(T) which are zero in a neighborhood of zeros(mo). 

The same argument can be applied to m' . But note that the right-hand side 
of l|5.2|) doesn't depend on m Q or m' Q (because v is the same). 

Therefore 



\m {z)\ 



rdis(z) 



m 

\m' (z)\ 



rdv(z) 



(5.3) 
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for all / 6 C(T) which are zero on a neighborhood of zeros(mo) U zeros(m ). 
From l|5.3|) it follows that \mo\ = \m' \, v- almost everywhere on 

C := T\(zeros(m ) U zeros(m )). 

Since the support of v is T, this implies that |mo| = on a set which is dense 
in T, and since the zeros of mo and m' are finite in number, |mo| = |m | on a 
dense subset of T. By continuity therefore 

|m | = \m' \ on T. 



If mo and W := \mo\ are not assumed continuous, there is still a variant of 
Theorem l5.ll but with a weaker conclusion. For functions W on T we introduce 
the following axioms. The a.e. conditions are taken with respect to the Haar 
measure /ionT: 

(i) W G L°°(T), 

(ii) W > a.e. on T with respect to fj,, 

(iii) i ^ tf/(w) = 1 a.e. on T, 

(iv) The limit 

n 

aV w := lim -Viy (fe) d/z (5-4) 

fe=l 

exists in Mi(T). Here, as before, 

W ik) (z) := W(z)W(z W ) ■ ■ ■ W (z^ 1 ) , (5.5) 

and 

(R w f)(z) = ± W(w)f(w), zeT, feC(T). (5.6) 

W N = Z 

Theorem 5.4. Let N G Z +7 and Ze£ a function W be given on T satisfying 
(i)-(iv) above. Let vw , W^ k ' and Rw be given as in (|5.4|) ~ (|5.6|) . Then (a)-(c) 
hold: 

(a) v w {Rw{f)) = vw{f), /GC(T). 
(6) If Wq, and W both satisfy (i)-(iv), set 

1 n 

di/ := lim - Viy (fe) aX- (5-7) 

fc=l 

/^^o (i?w(/Wo)) 

defines a measure on T which is absolutely continuous with respect to vq with 
Radon- Nikodym derivative W; i.e., we have 

vo(Rw{fW Q )) = MfW), feC(T). (5.8) 
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(c) If Wq and W are as in (b), then the following are equivalent: 
(cj) i/ % = v Q 
and 

(eg) There is a Borel subset E C T such that vq(E) — 1 and 

W (z) = W(z) (5.9) 

for all z € E . 

Proof. The structure of this proof is as that of Theorem l5.ll the essential step 
consists of the following two duality identities. Each one amounts to a basic 
property of the Haar measure /i on T. For every k £ Z + and / S C(T), we have 

Rw(f)W ik) dfi = J fW {k+1) d^ 

T T 



and 



/ Rw(fWo)W (k) d f i = J fWW^ k+1) d^. 



T T 

Using these, and (|5.2I) for the measure vq, the desired conclusions follow as 
before. For (c), in particular, we note that (ci) yields the identity Wdvo = 
Wodvo for the two functions Wo and W on T. Hence Wo and W must agree on 
a Borel subset in T of full i^o-measure, and conversely. | 

6 Transformation Rules 

If mo: T — *■ C is a Fourier polynomial; i.e., represented by a finite sum rao(z) = 
J2 k a kZ k , and if 

1 ]T |m H| 2 = l, zeT, (6.1) 

W N =Z 

for some N £ Z+, N > 2 we showed in |2j that there are functions mi , . . . , mjy_i 
on T such that 

jy m j{w)m k (w) = 8j,k, zeT, (6.2) 

W N =Z 

or equivalently, the N x N matrix 



(6.3) 



is unitary; i.e., defines a function from T into the group C/jv(C) of all unitary 
N x N matrices. Moreover the functions mi, . . . , mjsr-i may be chosen to be 
Fourier polynomials of the same total degree as mo- For the example in Section 
2, N = 3, and mo(z) — ^| , condition (|6.1|l is satisfied, and the other two 
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functions may be chosen as: mi (z) — z, and 7712(2) — 1 . There is a more 
general result, also from which defines a transitive action of the group Gat of 
all unitary matrix functions (Gat is often called the A^th order loop-group, and 
it is used in homotopy theory). An element in Gat is a function A: T — > Un(C). 
Let Fjy denote all the functions m = (mj)J^ 1 : T — > which satisfy l|6.2[) . or 
equivalently (|6.3() . Define the action of A on m as follows: 

m A (z) := A(z N )m(z), zeT (6.4) 

Lemma 6.1. The action of Gat on fjv is transitive and effective; specifically, 
for any two m and to' £ Fn, there is a unique A S Gjy such that to' = m ; 
i.e., A transforms to to m! . 

Proof. For functions / and j on T, 



1 



(f,g) N (z):=jj Yl /(«M«0, ^T. (6.5) 



„JV — , 



If to, to' G Fn are given, set 

Aj t k(z) := (m k ,mj) N (z), zel (6.6) 



Then an easy verification shows that A — {Aj.k) 1 ^ k l Q defines an element in Gat 
which transforms to to to'. Conversely, if to € Fn is given, and m A is defined 
by then it follows that m A <G F N if and only if A <G G 



AT- 



If a function too: T — * C is given, and satisfies (f Cj . 1 1> for some N, then the 
Ruelle transfer operator 

(i?/)(z):=l KMl'/M, (6.7) 

M> N =Z 

satisfies Rl = 1 where 1 denotes the constant function 1 on T. 

By a probability measure on T, we mean a (positive) Borel measure v on T 
such that f(T) =1. The probability measure will be denoted M\(T). 

Terminology for measures v on T: If / 6 G(T), set 

!/(/) := ffdv=f f{z)dv{z). 

T T 

If t: T — > T is a measurable transformation, set 

^ 1 (S) = Kr- 1 (^)) 

for Borel sets E cT. 

If i?: G(T) ->G(T) is linear, set 



{ V R){f) = v{Rf) = | (iZ/Jdi/. 
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If too is given as in l|2.13[) . we introduce 

£(m ) :={i/GMi(T) : vR^ 



v}. 



(6.8) 



If v G £(too) and £ C T is a Borel subset, we recall that E supports v 
if v(E) = 1. Note that from the examples in Section 4, it may be that the 
support of v is all of T even though v has supporting Borel sets E with zero 
Haar measure; i.e., v (E) — 1 and n(E) = 0. 

We now return to the case of the middle-third Cantor set C. Set s := log 3 (2), 
and view \c as an element in the Hilbert space L 2 (M, (dx) s ). We recall the usual 
unitary operators (Uf)(x) := -^/(f ), and (T k f)(x) — f(x — k),k£ Z, and the 
relation 

k £ Z. (6.9) 



unu- 1 = T 3k , 



If m (z) — a k z k , m £ L°°(T), is given, we define the cascade approxi- 
mation operator M = M mo as before 

Mf(x) = U^moiTjfix) = V2^a k f(3x - k), f £ L 2 (R, (dx) s ). 



The condition 



I l m oM! 2 = !j a - e - z£T, 



(6.10) 



will be a standing assumption on too. We then define the sequence rriQ (z) :— 
mo(z)mo(z 3 ) ■ ■ ■ mo(z 3 ), and we say that too has frequency localization if the 
limit of an associated sequence of measures, 



lim 



TOq (z) dfi(z) 



exists; i.e., if there is a v £ Mi(T) such that 



lim / f(z 



mf'(z) dfi(z) = / f(z)dv(z) 



(6.11) 



holds for all / £ C(T). Recall that if too £ Lip(T) is assumed, and R„ lQ has 
Perron-Frobenius spectrum, then it has frequency localization, and the limit 
measure v satisfies 



lim dist 



Haus 



[dv, 



» 



dfx) = 0. 



(6.12) 



Theorem 6.2. : The Dichotomy Theorem. Let too £ L°°(T) be given, and 
suppose ^(j.lUfi holds, and let v be the corresponding limit measure. Assume 
further that, for k £ 1+, 



lim / \mfr\z)\ 2 A^(z)d»(z) = v(A%>>) 



|2 



(6.13) 
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where A is the matrix function defined in j6'.6J) . and 

A^ k '(z) := A(z)A(z 3 )...A(z 3 ). Let M = M mo be the cascade approximation 
operator in L 2 (R, (dx) s ), s — log 3 (2). Then the limit 

lim M"xc exists in L 2 (R,(dx) s ) (6.14) 

n — >oo 

if and only if there is a Borel subset E C T such that v(E) — 1 (i.e., E is a 
supporting set for v), and Too (z) = , for all z 6 E. In the special case 
where A is further assumed continuous and too has frequency localization, the 
condition <6'.J51) is automatically satisfied, 

Mf(x) = (M c f)(x) = /(3s) + /(3s - 2) (6.15) 

and 

M cXc = Xc- (6-16) 

Proof. Suppose first that (|6.14|) holds; i.e., that the cascading limit exists in 
L 2 (R, (dx) s ). Then in particular 

lim ||M" Xc - M^ Xc \\ L2{{dxy) = 0. (6.17) 

We saw, using |Sj that there is a measurable matrix function A: T — > [/3(C) 
such that too is the first component in the product, matrix times vector, 



m (z) 

mi (z) ) =A(z 3 ) 
m 2 (z) 



\ 



z I . zeT. (6.18) 

1 — z 2 
V2 



If the functions Aj^ denote the entries in the matrix A on the right-hand side 
in (|6.18|l . we saw that 

P(Xc M Xc)(z) = A 0fi (z), z 6 T (6.19) 

and 

(M"Xc,^" +1 Xc) L2((d;c)s) = y R n {p{xc,M Xc ))dn = J R n (A^W (6-20) 

T T 

where 

R n (A , ) (z) = i ^ |TO n) H| 2 A ,oH, zeT. 

it; 3 ™ 

After a change of variable in (|6.13|) . we conclude from l|(j.20|l that 

Jim (M n x c ,M n+1 Xc ) L2{{dx)s) = f A 0fl {z)dv{z) =: u(A , ) (6.21) 
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and therefore 

= ^\\M n Xc-M n+ \ c \\l 2{{dx)s) =2-2Bep(Ao, ). (6.22) 

From (|6.18|) . we know that \Aq^\ < 1, pointwise for z E T. From this, we get 
that ^(^4o,o) = 1- Hence, there is a Borel subset £cT, such that v(E) = 1, 
and ^0,0(2) = 1 for all z G E. Since 

|^oW| 2 + |A a(z)| 2 + |^ 2 (z)| 2 = l, zeT, (6.23) 

we conclude that Aq^i(z) = Ao^(z) — for z £ E. Using (|6.18() again, we finally 
get mo(z) = for z G E; i.e., the conclusion of the theorem holds. If Ao.o is 
assumed continuous, then Aq^o — 1 on T since the support of v is all of T. The 
conclusions <|6.15|I - H6.16[1 in the theorem then follow. 

We now turn to the converse implication: If some supporting set E C T 
exists such that Aq$(z) = 1 for z G E, then j/(^4o,o) = / T Ao.odv = J E Ao^dv = 
j E dv — v(E) = 1. To prove the convergence in L 2 ((dx) s ) of the cascades in 
(|6.14|l . we must consider 

\\M n x C -M n+k Xc \\ 2 L2{(dx)s) = 2-2Re J R n (p( Xc , M k Xc )W, (6.24) 

T 

we note as in the case k = 1, that 

p(xc,M k Xc ) (z) = A$(z). 

If z G -E, then A ,o(^) = 1, and A 0j (z) = A,- )0 (.?) = for j = 1, 2. 
As a result, using l|6.23|l . we get 

2 

4 2 oW = E AwCO^oC* 8 ) = A ,o(z)Ao, (z 3 ) = 4>,o(A 
and therefore 

J A$(z)dv(z) = y A 0iQ {z 3 )dv(z) 

E E 

= J Ao,o{z 3 )dv(z) 

T 

= J A 0fl (z)dv(z) 

T 

= v(Ao,o) 
= 1. 

Continuing by induction, we find a supporting set, which is also denoted E, 
such that 

45(2) = 1 
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for all z £ E, k — 1, 2, • • • . Since 
substitution into l|{j.24[) yields 



\M n Xc ~ M n+k Xc 



,2 

\L 2 ((dx) 



2 - 2 Re / R 



T 

= 0. 



This proves convergence of the cascades, and concludes the proof of the theo- 
rem. | 

7 Low pass filters 

For functions on the real line R, and for every N E Z + , N > 2, the scaling 
identity takes the form 

<p(x) = a k ip(Nx - k) (7.1) 

fcez 

where D is a dimensional fixed constant. We take D — y/~N. The values a k are 
called masking coefficients, and 

m (z):=D- 1 J2akZ k (7.2) 

is the corresponding low-pass filter. The terminology is from graphics algorithms 
and signal processing, and the books [2] and [5] explain this connection in more 
detail. The function mo is viewed as a function on T = M./27rZ, or alternately 
as a 27r-periodic function on K, via z :— e~ 10 , 6 € R. It turns out that the 
regularity properties of too are significant for the spectral theoretic properties 
which hold for the operators associated with mo, specifically the cascade sub- 
division operator, and the Ruelle transfer operator. The function spaces which 
serve as repository for the function mo are measurable functions on T, for ex- 
ample L P (T), 1 < p < oo, the continuous functions; i.e., C(T), or the Lipschitz 
functions Lip(T). 

We will consider low-pass filters mo with the following properties: 

(1) to g Lip(T); 

(2) too has a finite number of zeros; 

(3) iZmoCl) = 1. 

Proposition 7.1. Suppose mo satisfies (1), (2), (3) and 

dim {ft G C(T) \R m „(h) = h} > 2. (7.3) 

Then there exists an (mo, N)- cycle (i.e., a set {zi,...z p } with zf — Zi+i, 
Zp = z\ and \mo(zi)\ 2 — N for all i). 
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Proof. Let h € C(T) satisfy R mo (h) = h, and ft non-constant. Taking the real 
or imaginary part, we may assume that h is real valued. Also, replacing h by 
Halloo — h 1 we may assume h > 0, and that /i has some zeros. 

We prove that all the zeros of h must be cyclic points. 

Suppose not, and let zq € T be a zero of h which is not on a cycle. Then, 
1 = zo and W2 2 = zo with l\ ^ £2 implies w\ ^ W2- Otherwise we have for 
some l\ < £2, Zq = w[ = zq so zq is a cyclic point. 

We say that w is at level £ if w N — z . The previous remark shows that £ 
is uniquely determined by w. 

Since h(zo) = 0, it follows that 

^ \m (w)\ 2 h(w) = 0, 

W N — Zq 

so I mo (w) 1 2 h(w) — for all w with uj^ = zo. Not all such w's can have 
mo(w) = because i? mo (l) = 1. Thus there is a z\ with zf = zo and ft(zi) = 0. 

By induction, there is a z n +i with z^ +1 = z„ and h(z n +i) = 0. 

Now, mo has only finitely many zeros, so from some level on, there are no 
zeros of m ; i.e., if hi" = z with I > £ then m (w) 7^ 0. But then look at 
those w's with w N = z^ . Since h(z£ ) — and mo(w) 7^ 0, it follows that 
h(w) = 0. By induction, h(w) = for all w with u; = zo and n large enough. 

However, these w's form a dense set because for any interval (a, i)cT there 
is an m big enough such that r m (a, 6) = T 9z£ (where r(z) = z w ) so there is 
a u; G (a,£>) with a) lv = zo. Since /i is continuous, it follows that h = 0, a 
contradiction. 

Thus, all zeros of /i are cyclic points. In particular zo, zi,...z„,--- are 
cyclic. 

Since z and zi are cyclic and z^ = z (hence they are on the same cycle), 
if w 7^ Zi, a;'' = z then iu is not cyclic so /i(w) 7^ 0, and therefore mo(w) = 0. 
But this implies (from i? mo (l) = 1) that |mo(zi)| 2 = N. We can do the same 
for all terms of the cycle generated by z\ and the conclusion of the proposition 
follows. I 

Remark 7.2. It is known generally that the dimension of the eigenspace in 
l|7.3|l depends on the metric properties of orbits in T under z — > z N where 
N > 2 is fixed. These orbits are called cycles. The solutions h to 

R mo (h) = h (7.4) 

are called |mo| 2 -harmonic functions and are important in the general theory of 
branching processes. Their significance for the present discussion is noted in 
There we prove that each solution h S i 1 (T), h > 0, h ^ 0, to equation l|7.4[l .- 
i.e., a non- negative harmonic function, is naturally associated with a system 
{H, U, T, f) where H is a Hilbert space, U and T are unitary operators in H 
satisfying UTU^ 1 = T N (i.e., equation (|4.1|) 1. and the vector ip e H, ip 7^ 0, 
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satisfies the general scaling identity 



Utp — mo(T)(p, (i.e., the abstract form of l|1.4|) or 114.14(1 .) 

See also QE]. 

T generates a representation of L°° (T) by 

tt(/) = /(T), (/ g L°°(T)). 

This representation satisfies the commutation relation 

Uni^U- 1 = n(f(z N j), {f G L°°(T)). 

Iterating the scaling identity one has 

U n ip = i:{m^ ) )ip, (n>0). 

Moreover, the system (iJ, {/, 7T, </?) is determined from /i in (|7.4|) up unitary 
equivalence and it is called the wavelet representation associated to (mo, h). 

Returning to the form p(-, •) in 1(4.2(1 . we note that h is related to the new 
data by the two formulas, 

h =p((p,cp) 

and 

J £(z N )p(U<p, Uip){z)dfi{z) = J £(z)h(z)dn(z), for all £ G C(T). 

T T 

The paper |2S] treats a general form of the problem l|7.4|l . and these au- 
thors state that the number of (mo , JV)-cycles on T equals the dimension of the 
eigenspace in 1(7.3(1 .' i.e., the space of continuous i? mo -harmonic functions. Re- 
call the points {zi} on an (m , iV)-cycle satisfy |mo(zi)| = V~N, Zi+i = zf , and 
Zk = zq, if k is the length of the cycle. The following example shows that this 
result from 23 is in need of a slight correction: Take for example mo(z) = 
N = 3, where there are no (mo, 3)-cycles. What is true is that, if the dimen- 
sion in H7.3fl is > l,or if there are eigenvalues A G T\{1}, then there must be 
(mo, iV)-cycles, and the arguments in |25| work: All the invariant measures are 
supported on cycles. But if the dimension in 1)7.3(1 is 1, then there may, or may 
not, be (mo, A r )-cycles. If not, then the invariant measures have support equal 
to T. If {1} is an (mo, iV)-cycle, then Dirac's Si is an invariant measure. 

We stress this distinction because it is central to explaining our dichotomy; 
i.e., explaining when a non-zero solution ip exists to the scaling equation ((4.14(1 . 
If the given filter mo has an (mo, iV)-cycle of length 1, then there is a scaling 
function cp in L 2 (R), and we are in the classical (non- fractal) case. If mo and 
N are fixed, but there is no (mo, AT) -cycle on T, then tp is instead in one of 
the 7i s -Hilbert spaces, < s < 1, with s depending on the chosen particular 
iterated function system (IFS). 
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Proposition 7.3. Let mo be a filter that satisfies (1), (2), (3). If there exists 
a A ^ 1 with \X\ = 1 and h £ C(T), h ^ 0, swc/i i/ia£ R mo ( n ) = A/i, i/ien i/iere 
exists an (mo, N)-cycle. 

Proof. We know that the peripheral eigenvalue spectrum of R mo is a finite union 
of cyclic subgroups of T (see section 4.5 in [H]). Hence A" = 1 for some n. 

Then R^ o (h) = h so R m (n)(h) = h and /i is not a constant. But m is 

Lipschitz, it has finitely many zeros and R («)(!) = 1 (the scale for mo is 
JV n ). 

Using proDOsition l7.il it follows that there is an (m Q n \ 7V™)-cycle; i.e., there 



exist points z\,...,z v on T with = Zj+i, z^ = z\, and m Q n \zi) 



But | mo | 2 < N (since i? mo (l) = 1) and this implies that 



|moO*i)| = |^o(2i v )|" = • ■ 
and therefore Zj will generate an (mo, N)- cycle. | 

Theorem 7.4. Suppose mo satisfies (1), (2), (3). Then exactly one of the 
following affirmations is true: 

(i) There exists an (mo, N)- cycle. In this case, the invariant measures v with 

v(R mo (f)) = v(f),(feC(T)) 

are atomic supported on the (mo, N)- cycles. The spectrum of R mo is computed 
in and JSf. The wavelet representation associated to (mo, 1) is a direct sum 
of cyclic amplifications of I? (R) ( see m$) ■ 

(ii) There are no (mo, N)- cycles. In this case there are no eigenvalues for 
Rmo \c(t) with |A| = 1 other then A = 1; 1 is a simple eigenvalue. There 
exists a unigue probability measure v on T which is invariant for R mo (i.e., 
v(R mo (f)) = vU) for f SC(T)). 

lim iC (/) = Kf) uniformly f G C(T) . (7.5) 

n — >oo 

Proof. When there are no (mo, 7V)-cycles, proposition 17.11 and 17.31 show that 
there are no peripheral eigenvalues other than 1 and 1 is a simple eigenvalue. 
The statements about v follow from theorem 3.4.4 and proposition 4.4.4 in jS] 
and their proofs. | 

Remark 7.5. When mo satisfies (1), (2), (3) and 1 is a simple eigenvalue for 
Rm Q |c(T) then the invariant measure v is unique and l|7.5|l holds (see [Hj and 
HQI). In the case of wavelet filters in L 2 (R), mo satisfies the extra condition 



m {z l 



= N 



m (l) = VN 

so {1} is an (mo, iV) -cycle. The measure v is simply the Dirac measure 5\. 
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Lemma 7.6. Let mo be a filter that satisfies (1), (2), (3). Assume in addition 
that 1 is a simple eigenvalue for R ma \c(i)- Consider the wavelet representation 
(H,U,ir, ip) associated to (mo, 1) as in remark \l .i\ 
Then for all £ G H with ||£|| = 1 and all f G C (T), 

lim {i | U- n n(J)ZTt) = v{f). 



Proof. First take £ of the form f = {/ m Tr(g)ip with m £ Z, g E C (T). Then, 
for n > m: 

(U- m n(g)<p | U- n Tr(f)U n U- m n(g)<p) 
tt (.g (>"- m ) m -^(z)) p | tt (/(z).g (^ m ) m "- m >(z)) p 



/(*) 



(n— m) 



iflWi a i4r ro) /(«)d/i. 



Since ||£| = 1, it follows that ||7r(<7)</s|| = 1 so J T \g(z)\ 2 d[i = 1. 
Also, from lim,^ i2J£ TO (/)(z) = !/(/) uniformly. 

Therefore 

lim <£ | U- n 7r(f)U n C) = / lff(^)| 2 K/)<^ = K/) 



Now take £ G H arbitrarily, ||£|| = 1. We can approximate £ by a sequence 
of the form mentioned before, with = 1. 

Fix e > 0. Then there is a j such that £j — £|| < (|) H/H^ and there is an 
n f such that, for n> n f . 



Then 



K^I^"7r(/)^%->-K/)|<3- 



|(£|C/-™7r(/)^ n 0-K/)| 
< |<e I U- n n{f)U n £) - <£ | u- n n(f)u n Q\ 

+ |(€ I U~ n ^f)U%) - I C/- n 7r(/)f/%.)| 
+ |(£.|C/-» 7r (/)C/™^|>-^(/)| 

<ll/ILlk-^l|HCll + II/IL 11^-11 IIC -^|| + |<e 



Theorem 7.7. Let mo,m 6e two filters satisfying (1), (2), (3) and suppose 
that I is a simple eigenvalue for R mo and R m ' on C(T). Let ^ and v 1 be the 
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invariant probability measures for R mo and R m ' respectively and let (H, U, it, ip), 
(H',U',Tr',(p') be the wavelet representations associated to (mo, 1) and (m , 1) 
respectively. 

If v 7^ v' then the two wavelet representations are disjoint. 

Proof. Suppose the representations are not disjoint, then there is a partial isom- 
etry W from H to H', where H and H' are the respective Hilbert spaces, and 
W j= 0. Take £ in the initial space of W, = 1; then ||Wf || = 1. Using lemma 
Owe have for all / e C (T) 

v{f) = lim (i | U- n n(f)U n 

n — >oo 

= lim (W£ | WU- n ir{f)U n £,) 

n — »oo 

= lim (W£ | U'- n ir'(f)U' n W£) 

n — ^oo 

= «/(/)• 

Thus v = v'. | 

Corollary 7.8. Let mo be a filter that satisfies (1), (2), (3), and suppose 1 is 
a simple eigenvalue for R mo on C (T) . Let v be the invariant measure for R mo 
and let {H, U, n, ip) be the wavelet representation associated to (too, 1). Suppose 
(p' € H is another orthogonal scaling function with filter to . Then 1 is a simple 
eigenvalue for R m i on C(T). If m' satisfies also (1), (2) then the invariant 
measure v' for i? m < is equal to the one for R mo , i.e., v' = v. 

Proof. Repeating the calculation given in the proof of lemma I7?fi1 we have 
v{f) = lim (U- m Tr(g)!p' | U- n 7T(f)U n (U- m 7T(9W)) 
= lim / \g{zf R n ^ m {f)dn 

T 

For all f,g€C (T), to € Z, J T \g\ 2 d[i = 1. Suppose h € C (T), h non constant 
with Rm> (h) = h. Then 

v (h) = J \g{z)\ 2 h(z)dn 

T 

for all g € C (T) with J T |(?(z)| 2 dfi = 1 then /i is constant. 
The last assertion follows directly from theorem 17. 71 | 

Corollary 7.9. Let m be a filter that satisfies (1), (2), (3) and suppose there 
are no (too, N)-cycles. Then the wavelet representation associated to (too, 1) is 
disjoint from the classical wavelet representation on L? (R) . 

Proof. Since 5i is not invariant for R mo , everything follows from theorem !7.7l | 
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